Driven classical diffusion with strong correlated disorder 
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We analyze one-dimensional motion of an overdamped classical particle in the presence of external 
disorder potential and an arbitrary driving force F. In thermodynamical limit the effective force- 
dependent mobility j-i(F) is self- averaging, although the required system size may be exponentially 
large for strong disorder. We calculate the mobility fi(F) exactly, generalizing the known results in 
linear response (weak driving force) and the perturbation theory in powers of the disorder amplitude. 
For a strong disorder potential with power-law correlations we identify a non-linear regime with a 
prominent power-law dependence of the logarithm of (J,(F) on the driving force. 

PACS numbers: 05.40. Jc,72.20.Ht,72.80.Ng 



We consider stationary diffusion in one dimension 
(ID), in the presence of a random disorder potential and 
a constant driving field. We show that the effective field- 
dependent mobility fi(F) is self-averaging, and calculate 
its dependence on the field F for various forms of corre- 
lated disorder. In particular, for a stationary diffusion in 
a strong disorder potential with power-law correlations 
at large distances [Eqs. (||), there is a wide in- 

termediate range of values of the driving force F where 
the logarithm of the effective mobility scales as a non- 
trivial power of F (Fig. |l|). With Coulomb- like (n — 1) or 
longer-range correlations, the effective mobility is a sin- 
gular function of F already at F — 0; the applicability re- 
gion of the linear transport is essentially absent (Fig. |^). 
These properties remain in the presence of weak interac- 
tion between the particles, introduced here at the level 
of the self-consistent Poisson equation which describes a 
Debye-like screening in the presence of strong disorder. 

Our results AHbU 7 to a number of systems where clas- 
sical diffusionEraoQ is the main mechanism of trans- 
port. Other applications include carrier diffusion in semi- 
conductor nanostructures, nano-scale thermoelectricsEm, 
transport in DNA0 and other biological systems. In par- 
ticular, the mean-field interaction model is directly appli- 
cable to ionic transport through cell membranesa, where 
one expects a number of parallel identical channels. 

Single-particle diffusion in ID is described by the 
Smoluchovsky equation 




FIG. 1: Logarithm of the effective mobility renormalization 
ln(/i(F)//i(0)) [Eqs. (^), (H)] with the correlation function 

g(x) = (1 + x 2 /£)- n/2 for strong disorder (V /T = 5) and 
large driving forces F. Dotted lines guide the eye with the 
slope of the intermediate asymptote ([14]) . Thin solid line is 
the analytic result ( |l5| ) for n = 1. 
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where U(x) is the external potential and f(t) is the ther- 
mal force with the correlator f(t)f(t') = 2Tr]o5(t — t'). 
The usual assumption is that the bare viscous friction 
coefficient r/o is determined by fast scattering events off 
phonons, short-range disorder, etc. The potential U{x) 
in Eq. (fn) is thus the part of the overall potential remain- 
ing after averaging over some distance scale; its precise 
value depends on the specific physical system. If we as- 
sume the carriers have charge e, we can also define the 
bare mobility in the absence of disorder, U(x) = —eEx: 




Ho = x/E = e/rjo. 



(2) 
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FIG. 2: As in Fig. [l]but for small driving forces F. Thin solid 
lines indicate the non-linear correction calculated analytically: 
linear in F for n > 1, proportional to F" for n < 1, and 
Eq. @ for n = 1. 
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The equation |l]) can be also rewritten as the transport 
equation for the average particle density n = n(x,t) and 
the particle current j = j(x, t), 

d t n + d x j = 0, j = -D d x n - r)Q 1 nd x U(x): (3) 

the diffusion constant Do is related to the viscous friction 
coefficient r/o by the Einstein relation, Dq = T/r)Q. 

The usual transport problem corresponds to station- 
ary diffusion in the presence of a random Gaussian po- 
tential V(x) and a constant driving force F, with the 
total potential energy U(x) = V(x) — Fx. In the case 
of a periodic potential V(x) = V(x + a), the stationary 
solutioru'ElB corresponds to a constant average current j 
with periodic boundary conditions, n(a) — n(Q). Then, 
by normalizing the density profile n(x) over the period, 
we can use the current to define the average drift veloc- 
ity, v = j /n = ja, as well as the effective viscous friction 
coefficient rj = F/v = F/ja, 



n 



aT 



dx 



x+a e [V(x')-V(x)+F(x-x')]/T 

dX ' i _ e -Fa/T ■ W 



Both integrations extend over the entire period, and, as 
it often happens in classical transport phenomena, for 
a sufficiently large a the effect of disorder becomes self- 
averaging [the required size can be large, see Eq. ( j2(j ) 
below] . In such cases we can replace Eq. (|4|) by its average 
over disorder. We assume that the disorder distribution 
is Gaussian with the correlations 

(V(x)) = 0, (V(x)V(x')) = V 2 g(x - x'), (5) 

where the local r.m.s. value Vq is taken as the measure of 
the disorder strength and the correlation function g{x) is 
defined so that g(0) = 1. Then, in the thermodynamical 
limit a — > oo, the effective mobility fi(F) = e/(r]}, and 



The disorder correlation function g(x) is expected to de- 
crease with x, remaining substantially different from zero 
over the distance of the order of the appropriate corre- 
lation length £. It is clear from the weak-disorder ex- 
pression (H) that for such finite-range disorder there is a 
distinct crossover force F% ~ T/£: while F < F$ have 
relatively little effect on the mobility, larger values of F 
begin to suppress the effect of disorder as large-scale po- 
tential valleys and hills gradually disappear. 

To analyze an analogous effect for strong finite-range 
disorder, note that for large Vq/T, A x = e~ 9 ^ v °l T is 
exponentially small for x < £. This effectively limits the 
integration in Eq. (0) to the region x > £, so that 



ln( M (F)//i(0)) ~ Fl/T, 



(9) 



where £ w ^ up to a logarithmic correction. Such a de- 
pendence on the applied field is analogous to the loga- 
rithmic susceptibilitytJ typical for systems with activated 
transport. Here it can be understood as the diffusion 
limited by far-spaced maxima of the potential, with the 
particles concentrated in the intermediate low minima; 
the applied driving force F effectively reduces the energy 
gap between the minima and the maxima and therefore 
has an exponential effect on the mobility. 

The expression (^) is valid qualitatively as long as the 
effect of the disorder remains large, fio/[i(F) ^> 1. The 
precise value of £ and the prefactor depends on the details 
of the disorder correlation function. For example, with 
the exponential correlation function g{x) = exp(— x/£), 
the integration (Q) can be done exactly in terms of the 
incomplete gamma function; the asymptotic form for 
V 2 /T 2 > max(l, X = F£/T) is 



KF) = (Vo/T) 2 * x>i 

m(o) r( x + i) 



(2tt X ) 



-1/2 e F«/T 



(10) 



^0 _ V 2 /T 2 £ 



dxe- Fx ' T e~^<l T 



In the small- F limit (but at the samepti 
gives the usual linear response resultcrli^ 



exp(y 2 /T 2 ), 



F > 0. (6) 
> T), this 

(7) 



which can be understood as the average of the activation 
exponent of the difference between the highest maximum 
and the lowest minimum of the potential. Generically, 
these would be in different parts of the sample and so 
the friction renormalization factors onto a product of the 
two averages {e v l T ){e~ v l T ) = exp(VQ 2 /T 2 ), independent 
of the form of the correlation function g{x). 

Similarly, the i stat ionary limit of the dynamical per- 
turbation theorytJta is restored by expanding Eq. (|6|) in 
powers of Vq/T and integrating the result by parts, 



Mo 



v 2 r°° 

l-^/ o dxe- F ^g'(x)+O(V i /T% 



(8) 



where f = f [l + ln(V 2 / (TF£))] . 

The disorder-induced transport non-linearity becomes 
even more pronounced for long-range potentials, e.g., 
those with power-law correlations at large distances. 
With long-range correlations the far-spaced maxima and 
minima of the potential are not entirely independent and, 
therefore, even a weak driving force may have a notice- 
able effect. Specifically, consider a correlation function 
with the asymptotic form 



(11) 



With strong enough disorder [g(x m i n )V 2 /T 2 > 1] the 
integral (^) will be determined by large x 7 in which case 
the expression can be rewritten approximately as 



« l n (a), l n (a)= [ dxe~ 

K F ) Jo 



(12) 



where a = (F^/T) n V 2 /T 2 can be small or large within 
the strong-disorder domain where Eq. (O) is applicable. 
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For sufficiently large a, the integration can be done us- 
ing the Gaussian approximation around the maximum at 

x = H 1 /^ 1 ', 



2-irxn \ 1/2 



1 



e -*o(l+l/n) a >X. 



(13) 



As a result, logarithm of /i(F) is proportional to a power, 



ln( M (F)M0)) ~ 



n/(n+l) 



(14) 



with a terge temperature-dependent length parameter 
ii = CM (V /T) 2 / n , C n = (n + l}J +1/n /n [cf. Eq. ©]. 
For small aCl, the integration ( |l2| ) can be done pertur- 
batively in powers of a a F n V a 2 for n < 1 or, using the 

identity T n (a) — J 1 / n {a 1 f n ), in powers of a 1 /" tx FV 2/ '™ 
for n > 1. For Coulomb disorder, n = 1, the result is 
expressed in terms of the Macdonald function, 



M (0)/ M (^)| n=1 nl^a) = 2a 1 ' 2 K 1 {2a 1 / 2 ). 



(15) 



For very small a<l, Ii{a) ~ 1 — a\n(e 1 ~ 2 ' 1 /a) , where 
7 w 0.577 is the Euler's constant; the correction is lin- 
ear in F up to a logarithm. Clearly, for strong Coulomb 
or longer-range disorder, n < 1, the mobility /i(-F') is a 
singular function of the driving force at F = 0; the linear- 
transport regime is essentially absent. These asymptotics 
are illustrated in Figs, [j], || for a model form of the dis- 
order correlation function g(x) = (1 + x 2 / £j[)~ n / 2 . 

Self-averaging. Our conclusions on the scaling of 
mobility in strongly-disordered diffusive ID systems are 
based on the average, Eq. (||). To analyze the sample- 
to-samplc fluctuations, consider the irreducible average 
((ri 2 )) = (i] 2 ) - {if) 2 of the effective friction 77 [Eq. (§)], 



((r? 2 )) 

v 2 



a T 2 J J 

lV (16) 



/ dz ( A x+y A x+aA '_„A 1 

' 2 Z 2 Z-\ TT— Z — 



where the correlator A z = e~ 9 ^ v ° ^ T is periodic under 
z — > z + a. Note that both A z and A J 1 enter Eq. (16). 
Therefore, unlike in the average (|J) , both short- and long- 
distance disorder correlations affect the variance of 77. 

For weak disorder, the expansion of Eq. (16) in powers 
of Vq/T <C 1 begins with the quartic term, 



«r/ 2 



VrfF 2 



[dxe- Fx / T [dye- F y/ T [dz 
Jo Jo Jo 

X 9z (2g z + 9z+x+y + gz+x-y - 2g z + x - 2gz+ y ) 

OT/4 pa pa 

= ±£L dzg, du e~ Fu / T {ug% u - g' z+u ). (17) 
al Jo Jo 

The integration is simplified in the limits of weak and 
large F; the combined result is 



((v 2 



2a 



ri poo 

^min(4&, T/F), £ 2 = J^dxg 2 (x). 



(18) 



Here, £2 is yet another correlation length, finite for short- 
range disorder and for long-range disorder with n > 1/2. 
Clearly, for weak disorder the variation of 77 is small and 
it is further reduced with increasing system size a. 

The situation is different for strong disorder, Vq 3> T, 
which causes an exponential renormalization of 77; large 
fluctuations are also expected. In this case A" 1 has a 
prominent maximum at the origin. Consequently, the in- 
tegral ( pf ) gets an exponentially large contribution from 
a vicinity of the point z — 0, x — y. Using the steepest 
descent method, we obtain 



«77 2 )) _ 2vr J F l2 e 4y o 2 /T 2 



Vo 



aV 2 



du ■ 



-2Fu/T e -4 g (u)V 2 /T 2 

g"(u)-g"(0) 



- (19) 



For not exceedingly large F the result is determined by 
values of u away from the origiiL Then, the denominator 
can be replaced by a constantly —g"(0) = 2/£ 2 , and the 
integral acquires precisely the form of Eq. <M). Generally, 



_ r ^ FT H 2V 2 /T 2 . 



2a V 2 



(20) 



while F in the prefactor can be small, it is assumed to be 
large on the scale of the system size, Fa/T 1. The nor- 
malization in Eq. ( |20| ) is chosen so that for short-range 
disorder C ~ 1 [cf. Eq. (j|)]. For power-law correlation 
tail, C = l n (2 n + 2 a)/[l n ( a )] 2 , with a = (F^/T) n V 2 /T 2 
[cf. Eqs. (|Tl]), (|T^)1. Overall, we conclude that the effec- 
tive mobility is a self-averaging quantity in the thermo- 
dynamical limit. Of course, the required system size can 
be large if the fluctuations are strong. 

We verified these conclusions by simulating diffusion 
in a ID short-range random potential (not shown). With 
periodic boundary conditions the viscous friction (^) 
could be obtained by averaging the time t it takes a parti- 
cle to travel over one period, Ax = a. As expected, with 
increasing a, the corresponding disorder average (t)/a 
approached the inverse of the average drift velocity. 

Stationary diffusion with weak interaction. The 
considered-problem differs from the canonical Kramers 
problemEIBB of over-the-barrier transport, where it is the 
dynamical equilibrium that establishes the exponentially 
different particle numbers in the "reservoirs" on the two 
sides of the barrier. Here, we consider a situation corre- 
sponding to a typical resistivity measurement in a macro- 
scopic sample where the total number of particles does 
not change with the applied field. Then, the macro- 
scopic current would be determined solely by the average 
drift velocity. It is important that the quantity is self- 
averaging, as the explicit disorder averaging would not 
be necessary for large enough samples. Previously, sta- 
tionary driven diffusion was cons idered , within the linear 
response with random disordemE^EJii 2 !, and al so Jar ar- 
bitrary F in deterministic periodic potentialsBaE3t3. 

In finite-size systems such a situation arises naturally, 
e.g., when diffusing particles are charged and the elec- 
troneutrality condition needs to be satisfied. The sim- 
plest case corresponds to the Debye mean-field screen- 
ing, where the potential in Eq. (H) is modified by the 
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self-consistent potential, U(x) — > U(x) +eip(x). Specifi- 
cally, we consider a ID Poisson equation, 

<p" = -47re(n - n), (21) 

as would be appropriate for diffusion in a 3D system 
with ID-modulation (layered disorder), a parallel bunch 
of identical DNA molecules, or electrostatically-coupled 
ionic cell channels. 

In the static equilibrium, F = 0, the coupled Eqs. ||), 
( pl| ) correspond to the non-linear screening problem; with 
weak disorder, Vq <C T, the Debye screening length 
is K , n 2 = Anne 2 /T. The linearized self-consistent 
screening problem can be also solved with the non-zero 
driving force, F > 0; the solution involves two screening 
parameters, s± = (k 2 + / 2 /4) 1/2 ±//2, where / = F/T w 
j/{Don) = X. Clearly, in the weak-interaction limit, 
/ 3> k, the shorter screening length si 1 k f- 1 = T/F 
is determined by the driving force, while the longer one 
diverges, sZ 1 ~ / /n 2 . 

With the driving force and a strong disorder, the prob- 
lem is forbiddingly complicated. However, if the interac- 
tion is weak, the additional potential would be small, and 
the screening equations can be linearized. To this end, 
it is convenient to eliminate the density n from Eqs. (g), 
( pl| ) and write the self-consistent equation for the scaled 
gradient of the screening potential e = eip'(x)/T, 

e" + e'(^-f)- K 2 s = r-ss\r^ K 2 (^ + X-f). 



This equation can be rewritten identically as 




For weak interaction, the last term in the r.h.s. is 
quadratic in small n 2 and can be ignored. The remain- 
ing equation should be solved for e with zero boundary 



conditions at infinity. The relation between the driving 
force F = JT and the diffusion current j = XnD is 
established from the self-consistency condition that the 
screening does not modify the net driving field. Equiv- 
alently, the disorder-averaged (e) = 0. Approximating 
s+ w /, after some algebra we again arrive at Eq. (jfy. 

Conclusions. We analyzed the stationary ID prob- 
lem of a driven diffusion in the presence of a random dis- 
order potential. For large systems and/or in the presence 
of an interaction fixing the number of particles trans- 
port should be described in terms of the effective mobil- 
ity [i. Strong disorder significantly reduces the mobility 
and leads to its non-trivial scaling as a function of the 
driving field F and the temperature. With finite-range 
disorder, the dependence pt(F) can be described in terms 
of the logarithmic susceptibility, Eq. (0) , a generic form 
for problems with activated transport. For a strongly- 
driven system with power-law disorder correlations the 
logarithm of the mobility scales as a power of the driv- 
ing force, Eq. (fL4|). The disorder effect is especially pro- 
nounced for Coulomb-like correlations [see Eq. (|l5|)]: the 
field-dependent correction to mobility is singular already 
at F = 0. 

The obtained explicit results are applicable for a num- 
ber of systems where diffusion is the main transport 
mechanism, and can be especially useful for character- 
ization of disorder distribution. Both quantum effects 
and particle-particle interaction would further modify the 
functional form of mobility. For the purposes of com- 
parison with experiment at low temperatures, it would 
be useful to obtain the corresponding results beyond the 
usual asymptotic limits of zero temperature and delta- 
correlated disorder. 
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